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FRACTAL WEYL LAWS IN DISCRETE MODELS OF CHAOTIC 

SCATTERING 

STEPHANE NONNENMACHER AND MACIEJ ZWORSKI 

o 

' Abstract. We analyze simple models of quantum chaotic scattering, namely quantized 

open baker's maps. We numerically compute the density of quantum resonances in the 
^ ■ semiclassical regime. This density satisfies a fractal Weyl law, where the exponent is gov- 

^ I erned by the (fractal) dimension of the set of trapped trajectories. This type of behaviour 

is also expected in the (physically more relevant) case of Hamiltonian chaotic scattering. 
Within a simplified model, we are able to rigorously prove this Weyl law, and compute 
quantities related to the "coherent transport" through the system, namely the conductance 
' and "shot noise" . The latter is close to the prediction of random matrix theory. 

> ' 

o 

^ : 1. Introduction 

, The study of resonances, or quasibound states, has a long tradition in theoretical, numer- 

^ I ical, and experimental chaotic scattering - see for instance jl] and references given there. 
CLi' In this paper we discuss the laws for the density of resonances at high energies, or in the 
r-| ■ semiclassical limit, and the closely related asymptotics of conductance, Fano factors, and 
i "shot noise" . Our models are based on a quantization of open baker maps [H [121 E] and 

we focus on fractal Weyl laws for the density of resonances. These laws have origins in the 
mathematical work on counting resonances [22] • 

'k^ \ In ^ we present the compact phase space models for chaotic scattering (open baker's 

■ maps) and their discrete quantizations. The numerical results on counting of quantum 

. , resonances showing an agreement with fractal Weyl laws are given in ^ In ^we discuss 
a model which is simpler on the quantum level but more complicated on the classical level 
(this model can also be interpreted as an alternative quantization of the original baker's 
map, see ^4.2p . In that case we can describe the distribution of resonances very precisely 
f ^4.2p . showing perfect agreement with the fractal Weyl law. We also find asymptotic 
expressions for the conductance and the Fano factor (or the "shot noise" factor). The 
fractal Weyl law appears naturally in these asymptotics and an interesting comparison 
with the random matrix theory is also made f ^4.3p . 

To put the fractal Weyl law in perspective we review the usual Weyl law for the density of 
states in the semiclassical limit. Let H{q,p) = p'^ + V{q) be a Hamiltonian with a confining 
potential V and let E he a nondegenerate energy level, 

(1.1) H{q,p)=E =^ VH{q,p)^0. 
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Assume further that the union of periodic orbits of the Hamilton flow on the surface H'^{E) 
has measure zero with respect to the Liouville measure. Then the spectrum of the quantized 
Hamiltonian, 

(1.2) H = -ffA + V{q) , qeW, 
near E satisfies, 

(1.3) i^^^Spec{H)n[E-ph,E + ph]} = j^ J 6{H{q,p) - E) dqdp + o{h-''^') , 

see jni for references to the mathematical literature on this subject. 

Suppose now that V{q) is nonconfining. The most extreme case is given by V{q) vanishing 
outside a compact set. An example of such a potential with g G is given in Fig. ^ In 
that case the eigenvalues are replaced by quantum resonances which can be defined as the 
poles of the meromorphic continuation of Green's function, G{z]q',q), from Im^; > to 
Im^; < 0. By Green's function we mean the integral kernel of the resolvent: 

(1.4) {z-H)-\{q')= ! G{z-q\q)u{q)dq, w G Cr(M") . 

We denote the set of resonances by Res(if). Near a nondegenerate energy level p.l|) we 
have the following bound (compare with ()1.3|) for a closed system) PJ: 

(1.5) # I Res(^) r}{[E-ph,E + ph] - z[0, ^h]) I < C(p, 7) . 

When the interaction region is separated from infinity by a barrier, this bound is optimal 
since resonances are well approximated by eigenvalues of a closed system. In that case the 
classical trapped set, 

(1.6) Ke = {{q,p) E H-\E) : $^(g,p) 00 , t ^ ±00} , 
has a non-empty interior in H~^{E), so that its dimension is equal to 2r2 — 1. 

Suppose now that the classical flow of the Hamiltonian H is hyperbolic on Ke, as is the 
case for instance in some energy range for the 2-D potential represented in Fig. ^ [[TH I19j. 
Following the original work of Sjostrand [12], the general upper bound ()1.5p is replaced by 
a bound involving the upper Minkowski dimension of Ke'- 

dimJ^^; = 2n — 1 — sup < c : limsup e""^ vol({p ET^e '■ d{KE,p) < e}) < 00 

We say that Ke is of pure dimension if the supremum is attained. For simplicity we assume 
that this is the case. Then under the assumption of hyperbolicity of the flow, we have [201 

(1.7) i^\^Res{H)n{[E-ph,E + ph]-i[0,-fh])^ <C{p,-f)h-''^ , 2pE + l = dimKE. 

This bound is expected to be optimal even though it is not clear what notion of dimension 
should be used for the lower bounds. The best chance lies in cases in which Ke has a 



FRACTAL WEYL LAWS IN CHAOTIC SCATTERING 



3 




Figure 1. A three bump potential exhibiting hyperbohc dynamics on a 
certain energy range. 

particularly nice structure. A class of Hamiltonians for which that happens is given by 
quotients of hyperbolic space by convex co-compact discrete groups . 

A fractal Weyl law for the density of resonances in larger regions are easier to verify and 
more likely to hold in general: 

(1.8) #I^Res{H)n{[E-6,E + 6]-i[0,-fh])^ ^C{6,-f)h-^^-\ 5 > fixed. 

The precise meaning of ~ is left vague in this conjectural statement. The exponent in 
this relation has been investigated numerically in a variety of settings and the results are 
encouraging 0. 

Acknowledgments. The first author thanks Marcos Saraceno for his insights and com- 
ments on the various types of quantum bakers. He is also grateful to UC Berkeley for 
the hospitality in April 2004. Generous support of both authors by the National Science 
Foundation under the grant DMS-0200732 is also gratefully acknowledged. 

2. Open baker maps and their quantizations 

We consider = [0,1) x [0,1), the two-torus, as our classical phase space. Classical 
observables are functions on and the classical dynamics is given in terms of an "open 
symplectic map" B, that is a map defined on a subset P C T^, which is invertible and 
canonical (area and orientation preserving) from V to B{T>). The points of \ P are 
interpreted as "falling in the hole" , or "escaping to infinity" 
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Following a construction performed in we will be concerned with open versions of 
the baker's map, obtained by restricting the "closed" baker's map to a subdomain of T^, 
union of vertical strips. As an example, if we restrict the 3-baker's map A^: 

( q' = 3q, p' = p/5, if0<g<l/3 
(2.1) A3(g,p) = (g',p')= < g' = 3g-l, p' = (p+l)/3, if 1/3 < g < 2/3 . 

[q' = 3q-2, p'={p + 2)/3, if 2/3 < g < 1 

to the domain 1^3 = \ {1/3 < g < 2/3}, we obtain the open 3-baker's map B^: 

This open map admits an inverse -83^^, which is a canonical map from B^iV-^) to V^. In 
this paper we will present numerical results for an open 5-baker's map, defined as 



(2.3) B,iq,p) = iq',p') 



def / g' = 5g-l, p' = {p + l)/5, if 1/5 < g < 2/5 
g' = 5g-3, p' = (p + 3)/5, if3/5<g<4/5. 



One can think of ^3 as model for a Poincare map for a 2-D closed Hamiltonian system. 
Removing the domain {1/3 < g < 2/3} from the torus corresponds to opening the system: 
the points in this domain will escape through the hole, that is, never come back to the 
Poincare section. In the context of mesoscopic quantum dots, such an opening is performed 
by connecting a lead to the dot, through which electrons are able to escape (see ^4.3|) . 

For open maps such as i? = i?3, we can define the incoming and outgoing tails, made of 
points which never escape in the forward (resp. backward) evolution: 

X G r_ C Vn > 0, 5"(x) G V . 

X G r+ ^ Vn > 0, 5-"(x) G B{V) . 

In the case of the map ()2.2|) . P^ = C3 x [0, 1), P+ = [0, 1) x C3, where C3 is the standard 
|-Cantor set on the interval (see Fig. |21). 

In analogy with p.fjj) . we also define the trapped set i^' = P+ fl P_ and, for any point 
X E K, its stable and unstable manifolds W±{x). In the case of the open 3-baker B3, we 
easily check that 

fi = dimP_ r]W+ = dimP+ n iy_ = - dim fsT = |^ . 

2 log 3 

The quantization of the open map ()2.2|1 is based on the quantization of the "closed" 
baker's map ^3. That, in an outline, is done as follows PJE]- To any iV G N we associate 
a space Hn — of quantum states on the torus. The components ipj, j G = 
{ 0, . . . , — 1 } of a state ip G Hn are the amplitudes of ip at the positions q = qj = 
(j + \)/^) and we will sometimes use Dirac's notation ipj = (gj|'0)- The choice of these 
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q 1 

Figure 2. We show, from left to right, approximations of the in/outcoming 
tails r_, r+ and the trapped set K for the open 3-baker B^. On the left and 
central plots, each color corresponds to points escaping at the same time. 



"half-integers positions" is justified by the parity symmetry g — > 1 — g they satisfy [THj, 
and is further explained in ^ The scalar product on Ti^r is the standard one on C^: 



(2.4) y^,(j)enN, (</)|^) = ^ 

j=0 

A classical observable depending on g G M/Z only, / = /(g), is obviously quantized as the 
multiplication operator 



The discrete Fourier transform 

(2.5) {gr,),,r = N~'/' e-2-(^-+^)(^-'+^)/^ , j, f = 0, . . . , N - 1 

transforms a "position vector" ipj = {qj\ip) into the corresponding "momentum vector" 
(PilV') = (^Af^)j- The momenta are also quantized to values = + /N , j = 0, . . . , N — 
1. Comparing the definition ()2.5|) with the (standard) Fourier transform on R, 



we see that the effective Planck's constant in the discrete model is ^ = (2ttN) ^. 
As a result, any observable g = g{p) can be quantized as 

Opjv(^7)^ = Q'k^MaiiJ + i/2)/iv)) Qn^. 
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The Weyl quantization on the torus generahzes this map / i— Op^(/) to any classical 
observable /, that is any (smooth) function on the torus, in such a way that a real observable 
/ is associated with self-adjoint operators, and 

^[Op^(/), Op^(^7)] = OpM,9}) + o{ff) . 

Let us now consider the following family of unitary operators on Ti-N, where is taken as 
a multiple of 3: 



(2.6) As,pos = As,N = G*N I g^/s 

Gn/s 

Since Qn exchanges position and momentum, the mixed momentum-position representation 
of is given by the matrix 

/ Gn/3 

^3,mom-pos = Gn ^3,Ar = I Qn/3 

V gN/3 

In terms of the quantized positions qj and momenta pk, the entries of this matrix are given 
by 

(^3,mom-pos)fc^. = (Pfcl^skj) = ^^|=^ exp ^-^i^Qj ' ^){Pk - g) J , 
I i+l i i+l 

- < < , 3 < Pfc < ' ^ = 0,1,2, 

and zero otherwise. One can then observe jl^ that for i = 0,1, 2, the function Se{p', q) = 
(3g — i){p' — i/3) generates the canonical map {q,p) ^— > {q',p') = (3g — i,p/3 + i/3) on the 
domain {q,p' G [i/3,{i+ l)/3)}, that is, the map ^3 ()2.1|) on this domain. The matrix 
elements {pk\As\qj) therefore exactly correspond to the Van Vleck semiclassical formula 
associated with the map A3. For this reason (and the unitarity of A3), the operator A3 was 
considered a good quantization of ^43 by Balazs and Voros. A more precise description of 
the correspondence between ^3 and ^3, including the role played by the discontinuities of 
A3, is explained in [HI §4.4]. 

To quantize the open baker B3 (j2.2p . we truncate the unitary operator A3 using the 
quantum projector on the domain V, U-z) == Op^{l.v) [Ej: in the position basis, we get 



Gn/s 



(2.7) S3,pos = B3,N = A3,N IIv = G*n\ I , AT e 3N . 

Gn/s 

This subunitary operator is a model for the quantization of a Poincare map of an open 
chaotic system jE]. The semiclassical regime corresponds to the limit N ^ 00. Similarly, 



Qn/3 
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the quantum open map associated with the 5-baker ()2.3|) is given by the sequence of 
matrices: 



^5,N — t'A 



N 



/ 











o\ 





Qn/5 



































Gn/5 





^0 











0^ 



G 5N. 



Let us now describe the correspondence between the resonances of a Schrodinger operator 
H, and the eigenvalues of our subunitary open quantum maps B^^n or B^^^ (denoted 
generically by B^)- 

Since -Bat is obtained by truncating the unitary propagator A^, it is natural to consider 
the family of truncated Schrodinger propagators X^~^^^^^Xi where x(q') is a cutoff function 
on some compact set supporting the scatterer. Although the precise eigenvalues of these 
propagators depend nontrivially on both x and the time t, these propagators admit a long- 
time expansion in terms of the resonances Zj of if [3]. At an informal level, one may 
write 

Zj€Res{H) 

On the other hand, the iterated open quantum map (B^)"' can obviously be expanded in 
terms of the eigenvalues Xj of Sat. For this reason, it makes sense to model the exponentials 



-izj/h 



by the eigenvalues Xj of our open quantum map B 



N- 



Upon this identification, the boxes in which we count resonances in (jl.7p . [E — ph, E + 
ph] — i[0,'yh], correspond to the regions 



(2.9) 



Aap = { 1 > |A| > r , I arg(Ae''')| < p } , r = exp(-7) G (0, 1) 



These analogies induce a conjectural fractal Weyl law for the quantum open bakers (j2.2l2.3j) 
which we now describe. 

First of all, we consider the partial dimension of the trapped set of the open map B: 

dimK 



p 



dim(r_ n 



Then, for any r G (0, 1), there should exist C{r) > (a priori, depending on the map B) 
such that, in the semi classical limit, the number of eigenvalues of Sat in the sectors ()2.9p 
behaves as 



(2.10) 



# { A G Spec(57v) n } 



271 



C(r) A^'^ A^ 



oo - 



The angular dependence p/ (27r) on the RHS means that the distribution of eigenvalues is 
expected to be asymptotically angular-symmetric. 
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In [T3j, the quantum 2-baker A2^n was decomposed into the block Q"^ 




:) 



and the 



complementary one. The spectral determinant for the unitary map, det{l—zA2^N), was then 
expanded in terms of these blocks. Although the classical open map associated with each 
block is quite simple (all points except a fixed one eventually escape) , the spectrum of each 
block was found to be rather complex, and quite different from semiclassical predictions. 

A scaling of the type ()2.10j) was conjectured in |TQ for another chaotic map, namely the 
open kicked rotator. This conjecture was then tested numerically, and a good agreement 
was observed. The scaling law A^^ was explained heuristically by counting the number of 
quantum states in an ^-neighbourhood of the incoming tail r_ (so that /i was effectively the 
dimension of r_). For the kicked rotator, the fractal exponent fi was not known analytically, 
and the authors related it to the mean dwell time of the dynamics, that is, the average 
time spent in the cavity before leaving it: ~ 1 — (Ardweii)"^ (A is the mean Lyapounov 
exponent). For our open baker's maps B^, the dwell time is Xdweii = 3, so the above formula 
is not valid: log 2/ log 3 7^ 1 — 1/(3 log 3). However, the above formula should give a good 
approximation of fi for a system with a large dwell time (that is, a small opening). 

In [jni we provide numerical evidence for the validity of ()2.10j) in the case of the open 
5-baker i?5 (j2.3p . at least when taking A^ along geometric subsequences. In ^4.21 we then 
construct a related quantum model, for which we can prove this Weyl law and calculate 
C{r) explicitly for inverse Planck's constants of the form N = 5'', k 



We numerically computed the spectra of several open baker's maps; in 12 , §5] we showed 
the numerical results concerning the 3-baker ()2.2|) . For a change, we will discuss here 
the open 5-baker ()2.3p . quantized in ()2.8|) . For this open map, the partial dimension of the 
repeller is /i = log2/log5 = 0.4306765... Compared to the 3-baker, this smaller exponent 
implies that the spectrum of B^ is expected to be much sparser than that of B^. For this 
reason, we will represent the spectra using a logarithmic scale (see Fig. E}, and consider 
regions Ar^-^^p for values of r ranging from r = 0.5 down to about r = 0.001. 

Let us now briefly explain the choice of "half-integer quantization" for the quantum 
positions and momenta ^21- The open map B^ is symmetric with respect to the parity 
transformation Il{q,p) = [l — q, 1— p): HoB^^ = B^oll. The choice of quantization is made 
so that the associated quantum map -B5 also possess this symmetry, that is, it commutes 
with the quantum parity operator 11 defined as Il\qj) = — 11 — qj) = —\qN-i-j)- We can 
then separately diagonalize the even and odd parts 



Both these operators have rank A^/5: together, they give the full nontrivial spectrum of 
Br,. We checked that the odd spectrum has the same characteristics as the even one, so we 



3. Numerical results 



,ev 



55o(l + n)/2 and 55,odd = 55o(l-n)/2. 
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only describe the properties of the latter. It is expected to satisfy the following fractal law 

(consequence of ()2.10|) ): 

(3.1) 

n{N,r) = #{Spec(55,^,ev) n A} ^ C(r) (iV/5)^°s2/iog5 ^ ^ ||^| > ^| ^ N ^ oo . 

The simplest set of iV's to test this fractal Weyl law is given by geometric sequences of the 
type {No X 5^ = 0, 1, . . .}: the law (jSlD means that the number of eigenvalues doubles 
when k ^ k + 1. In table Q we give some of the numbers n{N,r) along the sequence 
e {5^= X 20}, for some selected values of r. Along each column with r > 0.01, the 



= 20 X 5^= 


r = 0.5 


r = 0.1 


r = 0.05 


r = 0.01 


r = 0.05 


r = 0.001 


k = Q 


4 


10 


10 


13 


14 


16 


k = 1 


7 


19 


19 


25 


27 


35 


k = 2 


15 


36 


36 


48 


55 


122 


k = 3 


30 


69 


69 


104 


216 


402 



Table 1. Number of even-parity eigenvalues of B^^^ in {|^| > for N 
along the sequence {20 x 5*^}. 



numbers approximately double at each step k A; + 1, which seems to confirm the law 
(13.111 . The fact that this law fails for the small radii r = 0.005, 0.001 may be explained as 
follows: according to ()3.H1 . when A^ is large the huge majority of the A^/5 eigenvalues of 
-Bs AT^ev must be contained within an asymptotically small neighbourhood of the origin; if Ar 
intersects this neighbourhood, the law ()3.H) necessarily fails, since the counting function 
is proportional to A^ instead of A^'^. For the values of A^ listed in the table, this small 
region seems to be of radius > 0.005, explaining the departure from the fractal law in the 
last two columns. To further test the validity of the fractal law ()3.H) . we choose a set of 
values of r, and study the A^-dependence of n{N,r), for A^ taken along several geometric 
sequences, generalizing the above table. In Fig. |21 we plot this dependence in logarithmic 
scale for r = 0.3 (full lines), r = 0.1 (dashed lines), r = 0.015 (dot-dashed lines). Different 
geometric sequences are represented by a different colors. For almost all pairs {No, r) the 
points are almost aligned, and the slope is in very good agreement with the conjectured one 
H = log 2/ log 5. The less convincing data are the ones related to r = 0.3: for this radius, 
the numbers n{N,r) are still quite small, so that fluctuations are much more visible than 
for the smaller radii. We expect this effect to disappear for larger values of A^. 

In Fig. El the height of the curves does not only depend on r, but also on the sequence 
{No X 5'^} considered, especially for r = 0.3. To investigate this apparent contradiction with 
(13. ip we plot in Fig.Hlthe even spectra of i?5,Ar along 3 different geometric sequences. These 
plots suggest that, along a given geometric sequence, the eigenvalue density increases with 
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5 -baker along geometric sequences 

Even-parity states, logarithmic scale 



r=0.3 

-— r=0.1 




3 4 5 6 7 8 9 

log(N/5) 

Figure 3. Checking the A^-dependence of n{N,r) for various values of r 
and N along geometric sequences {No x S'^}. We also show (thick line) the 
conjectured slope log 2/ log 5, and give for comparison the slope 1 (thick 
dashed) . 

N uniformly with respect to (j) = argA, but very nonuniformly with respect to |A|. We see 
that some regions {ro < |A| < ri} remain empty even for large values of A^. The presence 
of gaps was already noticeable when comparing the second and third columns of table ^ 
obviously, for = 20 x 5^, there were no eigenvalues in the annulus {0.05 < |A| < 0.1}, 
which is confirmed visually in Fig. |3] (bottom). This non- uniform dependence on |A| implies 
that the profile function C{r) is nontrivial. 

The spectra for the two other geometric sequences also shows the presence of gaps, but 
the gaps differ from one geometric sequence to the other. This observation also contradicts 
the law 1)3.11) . In spite of these problems we nevertheless attempt to compute the profile 
function C(r) appearing in ()3.1|) . Fig. El (left) shows n{N, r) as functions of r G (0, 1), for 
along the same three geometric sequences (each one corresponding to a given color/width). 
We then rescale the vertical coordinate of each curve by the factor (A^/5)-i°s2/iog5^ ^^^^ 
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Figure 4. Even-parity spectra of the quantum baker's maps B^^n, along 
geometric sequences N — N^x 5'^ . The eigenvalues are represented on a loga- 
rithmic scale (arg A/27r against log |A|). The dimensions indicated correspond 
to N/5. 
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Figure 5. On the left, we plot n{N,r) as functions of r G (0,1). The 
numbers in the legend are On the right we have multiplied n{N,r) by 

the factor (iV/5)-'°g2/i°g5. 



plot the rescaled curves in Fig. El (right). From far away, these rescaled curves are fairly 
superposed on each other, which shows that the conjectured scaling (|3.1|) is approximately 
correct. Yet, a closer inspection shows that a much better convergence to a single function 
occurs along each individual geometric sequence. For instance, the curves for = 8 x 5'^ 
"pointwise" converge to the last one along this sequence (A^/5 = 5000), which has a plateau 
on {0.2 < r < 0.4} corresponding to a spectral gap. The curves of the two other sequences 
seem to converge as well, with plateaux on different intervals. 

In the case of the open kicked rotator studied in ^| the rescaled curves n{N, r) are 
more or less superposed, therefore defining a profile function C (r) . The authors claim that 
this function corresponds reasonably well with a prediction of random matrix theory |24j . 
Our results for the 5-baker's map contradict this universality: there does not seem to be a 
global profile function C(r), but a family of such functions, which depend on the geometric 
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sequence {No x S'^}, which could be denoted by C{No, r). The law ()3.1|) needs to be adapted 
by restricting to geometric sequences, which yields the following empirical scaling law: 

For any A^o > and r G (0, 1), there exists C{No, r) > such that, for N along {No x S'^} 
and k ^ oo, 

(3.2) # { A e Spec(SAr) n A^p } - IT C{No, r) ^ , 

Zn 

where ^r-,i?,p is given by ()2.9|) . In Fig.|31the different profile functions are uniformly bounded, 
Ci{r) < C{No,r) < C2{r), for some envelope functions < Ci(r) < C2(r). 

This weakening of ()2.10j) to geometric sequences makes sense for baker's maps of the form 
i?3, i?5, which each have a uniform integer expansion factor, leading to number-theoretic 
properties. In the case of a nonlinear open chaotic map (as the open kicked rotator of [iSj), 
there is no reason for geometric sequences to play any role, so we expect ()2.1()j) to hold in 
that case. 

4. A COMPUTABLE MODEL 

Because we are unable to analyze the spectra of the quantum bakers rigorously, we 
introduce simplified models. In the case of the 3-baker, we observe (see Fig. IHl left) that 
the largest matrix elements are maximal along the "tilted diagonals" 

(4.1) (n,m) = (3/+e,/+£A^/3), with / G {0, . . . , Ar/3-1}, £e{0,2}, eG{0,l,2}. 

These "diagonals" correspond to a discretization of the the map B^ projected on the po- 
sition axis. Away from them, the coefficients do not decrease very fast due to the Gibbs 
phenomenon (diffraction). The elements on the "diagonals" have moduli l/-\/3 + 0{1/N) 
and their phases only depend on i, e in the above parametrization. Our simplified model 
is obtained by keeping only the elements on the "diagonals" (see Fig. IHl right), set their 
moduli to 1/a/3 and shift their phases by 7r/2 (for convenience). Using the parametrization 
(jHH), we get: 

~ 1 2?7r 

(4.2) iBs,N)nm = -j^ exp (--(e+ l/2)(£+ 1/2)) . 

\ 

uj^l^ 

1 





1 



u^/^ 

1/ 



For = 9 and using u = e^'^*/^, the matrix reads 



B 



UJ 



1/4 



3,9 



V3 



/ 1 


































UJ 




















1 

















0.1/2 

















UJ 




















1 

















0.1/2 











^0 





UJ 
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The matrix B^ ^ can obviously not be considered as a "small perturbation" of B^^^, since 

3 — baker (S) 3 — baker (W) 



0.5 m/N 1 0.5 m/N 1 




N= 27 N= 27 

Figure 6. Matrices -83,27 (left) and its toy model -83,27 (right). The gray 
scale represents the modulus of the matrix elements (white= 0, black= 1). 



we removed many nonnegligible "off-diagonal" elements. Actually, by acting with B^^^ on 
Gaussian coherent states, one realizes that these matrices do not quantize the open 3-baker 
B3 fl2.2|l . but rather a more complicated multivalued map B3, built upon B3 as follows: 

1 

(4.3) V(g,p)GP3, B,{q,p)= \J{B,{q,p) + {0,j/3)}. 

i=-i 

We refer to jT2l Proposition 6.1] for a precise statement. As opposed to -B3, the multivalued 
map ()4.3j) is no longer obtained by truncating a canonical transformation, but it comes 
from three different transformations. -B3 can be considered as a model of propagation 
with ray splitting. Another interpretation is given by considering a Markov process with 
probabilities P{x',x) being allocated at each step to the image points x' of x. Explicitly, 
the probabilities take the form 

so that for each a; G P3, the sum of the weights associated with the three images of x is 
indeed 1. 

Some of the characteristics of the dynamics remain the same as for B3. The local dy- 
namics of each branch is the same, and -B3 sends all points in \ P3 to infinity. One can 
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define incoming and outgoing tails for B^, (see As opposed to tlie case of B^, tliese 
tails are not symmetrical any more: r_ = C3 x [0, 1), r_|_ = T^. Yet, these formulas are 
slightly misleading. The second one comes from the fact that any point x G has two 
preimages through B3, namely xo = (5'/3, 3p), X2 = ((g + 2)/3, 3p), so no point ever escapes 
to infinity in the past. However, to these preimages are associated the respective weights 
P{x,Xo), P(x, X2), the sum of which is generally < 1: there is thus a loss of probability 
through B^^, which is not accounted for by the definition of r+. 

In the next section we will show that the matrices -Bs^at can nonetheless be interpreted 
as quantizations of the original open baker B^, as long as one switches to a different notion 
of quantization, derived from a different type of Fourier transform (the Walsh-Fourier 
transform) . 

Families of unitary matrices A2^n with a structure similar to -Bs.at have already been 
proposed as an alternative quantization of the 2-baker's map A2 PHI- These matrices can 
also closely related with the "semiquantum bakers" introduced in In the context of 

quantum graphs (a recently popular model for quantum chaos), unitary matrices similar 
with A2,N (but with random phases) occur as "unitary transfer matrices" associated with 
binary graphs [21] • In this framework, the matrix A2^n would correspond to a graph with 

very degenerate bond lengths. In this framework, the matrix B^ i^ is directly related with a 

I ^ 1 2 

classical transfer matrix defined by = | (-S3,Ar)jj'| , which represents the classical 

Markov process on the graph. In our case, this transfer matrix is the discretized version of 
the transfer (Perron- Frobenius) operator associated with the open map B^. 

4.1. The Walsh model interpretation of B^^^. In this section, we represent the ma- 
trices B3 N in a way suitable for their spectral analysis. This can be done only in the 
case where iV is a power of 3. This representation is connected with the Walsh model of 
harmonic analysis. 

The latter originally appeared in the context of fast signal processing [8 . The major 
advantage of Walsh harmonic analysis (compared with the usual Fourier analysis) is the 
possibility to strictly localize wave packets simultaneously in position and in momentum. 
For our problem, this has the effect to remove the diffraction problems due to the discon- 
tinuities of the classical map, which spoil the usual semiclassics [T3]. 

A recent preprint fTOj analyzes some special eigenstates of the "standard" quantum 2- 
baker, using the Walsh-Hadamard transform (which slightly differs from the Walsh trans- 
form we give below) as a "filter". We are doing something different here by constructing 
our simplified model B^^n from the Walsh transform, as -Bs^at was constructed from the 
discrete Fourier transform (see 



^M. Saraceno, private communication. 
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We first select the expanding coefficient of the baker's map, which we denote by G N 
(the map ()2.2|) is associated with D = 3, the map ()2.3|) with D = 5). Once this is done, we 
will restrict ourselves to the values of along the geometric sequence {N = , G N}. 
In this case, the Hilbert space can be naturally decomposed as a tensor product of k spaces 
C^: 

(4.4) Un = {C'')i®{C%®---®{C%. 

This decomposition appears naturally in the context of quantum computation, where each 
(j^D represents a "quantum D-git", that is, a quantum system with D levels. Here, we 
reahze this decomposition using the basis of position eigenstates \qj) of T-Cn (see [THj for the 
case D = 2). Indeed, each quantum position qj = {j + l/2)/N, j E Z^k = {0, . . . , — 1} 
is in one-to-one correspondence with a word e = eie2 ■ ■ ■ made of symbols (D-gits) ee = 
Q(i) e Zd: 

k 

(4.5) j = J2eeD'-'. 

e=i 

The usual order for j G Zjyk corresponds to the lexicographic order for the symbolic words 
e G (Z^))'^. Associating to each -D-git a D-dimensional vector space (C^)^ with canonical 
basis {eo, ei, . . . , en-i}, the position eigenstate Ig^) G Hn can be decomposed as 

(4.6) \qj) = ® (S) ■ ■ ■ ® e^^, . 

This identification realizes the tensor product decomposition ()4.4|) . 

The Fourier transforms (j2.5|) . and the simpler one without the 1/2 shift, 

„-2injj'/N 

(4.7) {J^N)n' = y=r- , J,J' eZ^, N = D\ 



'N 

are defined by applying the exponential function x i— >■ e~^*'^^ to the products 
. k 

^=5^ D-^im{jf), where e^(jj') = q(j>Hj') • 

m=2-k i+i'-m+k 

If we replace in fl4.7|) the exponential e~^*'^^ by the piecewise constant function ej:,{x) = 
exp(— 2i7r[Z}x]/-D), and replace each em{jj') by its value em{jj') modulo D, we obtain the 
matrix element 

k 

(V.),,. D-^'^ en{ D~-en.Uf)) = D-''^ ^xp ( - ^ e,Uf)) 
(4.8) ^ ""^^'^ 

2m 



= n^~^/^exp(-^^Q(j)e,+i_,(/) 

£=1 

The matrix Vh defines the Walsh transform in dimension D^. 
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Because we have used the "half-integer" Fourier transform ()2.5p to define our quantum 
baker's map, we will need a slightly different version of Walsh transform, namely 

k 

im^' = llD-'/' exp ( - ^ + 1/2) (6,+i_,(j') + 1/2)) . 

i=i 

Both Vfc and Wk preserve the tensor product structure ()4.4|) : for any Vi, . . . ,Vn G C'^, 

(4.9) Vkivi ^---^Vk) = J^DVk ® ■ ■ -^J^dVi , >Vfc(fi 'S)---0Vk) = GnVk ® ■ ■ -^GdVi ■ 

These expressions show that Vk and Wk are unitary. 

Specializing the computations to D = 3, we are now in position to define the toy model 
B^^N (in the case = 3^^) as the "Walsh-quantization" of the 3-baker ()2.2p (as opposed to 
the "standard" quantization of the multivalued map 53 ()4.3j) ). Indeed, one can check that 
the matrix ()4.2|) can be expressed as 

f Wk-i 

(4.10) B3N = W*k\ 

V Wk-i 

This formula is clearly the Walsh analogue of the definition ()2.7|) of the "standard" quantum 
open baker B^ j^. From this definition and (j4.9|) . we see the action of -Bs^at on tensor 
products: 

(4.11) Bs,Nivi ■ ■ ■ ®Vk) = f2 ® ■ ■ ■ ® ffc ® G3'^0,2Vl , Vj e , 

where 7ro,2 is the orthogonal projector (in C^) on Ccq © Ce2. 

4.2. Distribution of resonances. Using ()4.11|) we can explicitly describe the spectrum 
of B^^N for = 3*^. The computation is identical with Section 6.2], so we only give 
the results. The generalized kernel of B^ j^f is spanned by the position states \qj) such that 
ee{j) = 1 for at least one index 1 < i < k. This corresponds to positions qj "far" from the 
Cantor set C3, so that the classical points {qj,p) are sent to infinity at a time n < k. This 
kernel has dimension 3^= - 2*^ = - A^i°g2/iog3_ 

The nonzero eigenvalues of B^^n are given by the set (see Fig. [7j) 

{A+} U {A_} U U |J{e2-^/'= A^T^/' A^/'} , where A+ = 1, A_ = ^ . 

i=o p=i V 3 

For each p E {1, . . . , /c — 1}, the k eigenvalues of modulus |A_|^/'^ = 3-p/2'^' asymptotically 
have the same degeneracy (^^/k as k ^ 00 (semiclassical limit), which shows that their 
distribution is circular-symmetric. Taking these multiplicities into account, we obtain the 
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0.4 I ' 1 ' 1 1 1 

-0.4 -0.2 0.2 0.4 

Figure 7. Eigenvalues of the matrices B^^n for = 3^° (triangles) and for 
N = 3^^ (crosses), forming lattices in a logarithmic scale. The two horizontal 
lines correspond to the spectral radius \z\ = 1 (thin line) and the "peak 
multiplicity" \z\ = 3~^/^ (thick line). Notice the difference of vertical scale 
compared with the spectra of Fig. |3J 

following Weyl law for the eigenvalues of B^^n inside a region ()2.9|) . along the sequence 
N e {3''}, k oo: 

# I Spec(53,iv) n Ar,s,p 1=7^^^ r) + o(l)) 

(4-12) log 2 

= dim(r„ n iy+) = ^ , C(l, r) = 1(0,3-1/4, (r) . 

The values A_, A+ in ()4.12j) are the nonzero eigenvalues of the matrix o 7ro,2 appearing 
in ()4.11|) . We used the notation C(l,r) for the profile function to be consistent with our 
notations in ()3.2p . that is, to emphasize that this estimate is valid only along the sequence 
G { 1 X 3^= }. 

We notice that the spectrum of the classical transfer matrix Bs^n defined at the end of 
12] is drastically different: this matrix admits one simple nontrivial eigenvalue A = 2/3 
(interpreted as the classical escape rate), the rest of the spectrum lying in the generalized 
kernel. Therefore, the features of the quantum spectrum is intimately related with the 
oscillatory phases of B^^n along the "diagonals" . 

4.3. Conductance and Shot Noise. In this section, we consider an open baker's map as 
a model of quantum transport through a "chaotic quantum dot" , that is a 2-dimensional 
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cavity connected to the outside world through a certain number of "leads" carrying the 
current; each lead is connected to the cavity along a segment Lj of the boundary (see 
Fig. IHl), and the connection is assumed to be "perfect": a particule inside the cavity which 
hits the boundary along q G Lj is completely evacuated to the lead. Therefore, the phase 
space domain Lj x [0, 1) above this segment is a part of the "hole", in the terminology of 
^ whereas the remaining set I = [0, 1) \ (UL^) represents the boundary of the quantum 
dot, which lifts to the phase space domain V = I x [0, 1). 

In the previous sections we have studied the open quantum map obtained by projecting a 
unitary quantum dynamics (called generically Un) onto a subdomain V of the phase space: 
resonances were defined as the eigenvalues of Un^v- These resonances are supposed to 
represent the metastable quantum states inside the open quantum dot, after it has been 
opened. In the present section, we want to study another aspect of the open system, namely 
the "transport" through the dot, using the formalism of [22] • We will focus on the case 
where the opening L splits into two segments L = LiU L2, and we study the transmission 
matrix from the lead Li to the lead L2. Once we are given, on one side, the quantum 




Figure 8. Closed and open chaotic cavity. On the left, a (bounded) trajec- 
tory is schematically associated with a power of the unitary quantum prop- 
gator (n represents the number of bounces). On the right, a "transmitting" 
trajectory is associated with a term of the matrix ()4.13|) . 

map Uj\f associated with the closed dynamics inside the "cavity", on the other side, the 
projectors on the leads 11^. and on the "interior" 11/ = U-jj, the transmission matrix (from 
Ll to L2) is defined as the block 

(4.13) m = ^e^'^'^n,., Um {UjUNr-'UL, . 

n>l 

The parameter G [0, 27r) is the "quasi-energy" of the particles. According to Landauer's 
theory of coherent transport, each eigenvalue Tj('(9) of the matrix t{-^)t*{-t}) corresponds to 
a "transmission channel" . The dimensionless conductance of the system is then given by 



(4.14) 



g{^) =tr {t{^)t*i^)) 
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A transmission channel is "classical" if the eigenvalue Tj is very close to unity (perfect 
transmission) or close to zero (perfect reflection). The intermediate values characterize 
"nonclassical channels" (governed by strong interference effects). The number number of 
the latter can be estimated by the noise power 

(4.15) P{^) = tr (^t{^)r{^) [Id - t{^)t*{^))^ , 

or equivalently the Fano factor, F = P/g. It is sometimes necessary to perform an ensemble 
averaging over to obtain significant results [22j. However, for the model we study here, 
these quantities will depend very little on {}. The closed quantum dot will be modeled by 



L 



L 



Figure 9. 4-baker's map modeling the chaotic cavity. The left- and right- 
most vertical strips correspond to the two openings ( "leads" ) . 



the following quantum map: we consider the 4-baker's map A4 and quantize it using the 
Walsh transform Vk ()4.8|) with D = 4. In dimension = 4^, our unitary propagator is 
therefore 














\ 





Vk-i 
















Vk- 


-1 














1 / 



We attach the leads on the intervals Li = [0, 1/4] and L2 = [3/4, 1]: this way, the projectors 
Uii as well as the projector 11/ = Id — Ul^ —11^2 can be represented as tensor product 
operators: 

Ilii = TTo ® Id4 (g) Id4 ® ■ ■ ■ , = VTs (g) Id4 (8 ■ • ■ , II/ = TT/ (g) Id4 (g ■ ■ ■ . 

Here tTj is the orthogonal projector on the basis state of C'^, and vt/ = tti © 712- This 
tensor action, together with the action of A4^]\f (analogous to (|4.1ip ). allow us to compute 
all quantities exphcitly. 

The spectrum of the "inside" propagator for this model, B4 N = jyH/, satisfies a fractal 
Weyl law of the type (j4.12j) along the sequence = 4^, with exponent fi = log 2/ log 4 = 
1/2, and profile C(l,r) = ]1[q 2-3/4] (r). 
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For this model and in the semi classical limit A; oo, we could compute the dimensionless 
conductance ()4.14|) . The computation [T21 §7.2] requires to control the time evolution up 
to n = Ck for some 1 < C < 2: this is of the order of the Ehrenfest time Te = k for the 
system. For any d we obtain 

(4.16) g{N = 4^ ^) = ^(1 + 0(2-'=)) = ^ (1 + o(l)) , A: - oo . 

Here N/A is the number of transmission channels from Li to L2, that is the rank of the 
matrix t{;d). We see that, as could be expected, approximately one half of the scattering 
channels get transmitted from one lead to the other, the other half being reflected back. 

Asymptotics for the shot noise (j4.15|) (which counts the "nonclassical" transmission chan- 
nels" ) are more interesting and again independent of d: 

(4.17) P{N = 4^ ^) = 2'-' Q + 0{e'^')) = 11 (iV/4)^ (1 + o(l)) , k 00 . 

Here fi = 1/2 is the dimension appearing in the fractal Weyl law for the resonances. A 
similar fractal law for the shot noise had been observed in [221 ^^^^ quantum 

kicked rotator; the power law A^^ for the number of nonclassical channels was explained 
there through a study of the dynamics up to the Ehrenfest time. 

The constant 11/80 in ()4.17|1 gives the average "shot noise" per nonclassical transmission 
channel. This number is close to the random matrix theory prediction for this quantity, 
namely 1/8 [71I22]- The precise number 11/80 certainly depends on which baker's map one 
starts from, and which quantization one uses. For instance, we did not check whether the 
"half-integer" Walsh quantization of the 4-baker leads to the same prefactor, but we expect 
the result to be close to it. It would be interesting to actually check the full distribution 
for the transmission eigenvalues Tj, and compare it with the prediction of random matrix 
theory [7j. 

The near agreement with random matrix theory is in contrast with the fact that the 
semiclassical resonance spectrum of the propagator B^^n inside the dot is very different 
from that of a random subunitary matrix. Somehow, the matrix t{'d), obtained by summing 
iterates of B^^n, has acquired some "randomness", as far as the distribution of its singular 
values is concerned. 

The transport properties of chaotic cavities has also been studied within the framework 
of quantum graphs. The shot noise (j4.15|) could be semiclassically estimated in the case of 
a "star graph", by summing over transmitting trajectories on the graph ^7] (they studied 
the case of "small openings"). The authors show that one needs to take into account 
the "action correlations" between different trajectories, in order to reproduce the random 
matrix result. As mentioned before, the matrix A4^i^ can be interpreted as the unitary 
transfer matrix for a different type of graph with bonds having degenerate lengths. 
Somehow, our use of the tensor product structure implicitly takes into account the action 
correlations for this particular graph. 
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5. Conclusions 

Quantum open baker's maps provide a simple and elegant model for the study of quantum 
resonances of open chaotic systems. The numerical investigation of these models is easily 
accessible and, as shown in ^ gives a good agreement with the fractal Weyl law on "small 
energy scales" , which is ()1.7|1 in the case of Hamiltonian flows. Only larger energy scales 
()1.8|) were considered previously. It would be interesting to investigate the spectrum of the 
model operator ()2.7|) for higher values of N . The naive numerical approach we took 

(full diagonalization of the matrices B^) only allowed to reach values < 5000. It would 
make more sense to use an algorithm allowing us to extract only the largest eigenvalues 
(which are the ones we are interested in), instead of the full spectrum. 

By modifying the standard quantum baker's map, in a way which still fits in the frame- 
work of quantization of chaotic dynamics, we obtained a model for which the fractal Weyl 
law ()4.12j) can be rigorously proven. Since the spectrum of this model is explicitely com- 
putable (and forms a lattice), it is forcibly nongeneric. However, the explicit computation 
of other physical quantities associated with our model, namely the conductance and the 
"shot noise", shows more generic properties. The fractal Weyl law is also present in the 
calculation of the "shot noise" , and the prefactor is (unexpectedly) close to random matrix 
predictions. 
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